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Abstract: Kasteleyn counted the number of domino tilings of a rectangle 
by considering a mutation of the adjacency matrix: a Kasteleyn matrix K. In 
this paper we present a generalization of Kasteleyn matrices and a combinatorial 
interpretation for the coefficients of the characteristic polynomial of KK* (which 
we call the singular polynomial) , where K is a generalized Kasteleyn matrix for 
a planar bipartite graph. We also present a (/-version of these ideas and a few 
results concerning tilings of special regions such as rectangles. 

Introduction 

Kasteleyn ([K]) counted the number of domino tilings of a rectangle by considering a 
mutation of the adjacency matrix, since then known as a Kasteleyn matrix ([LL], [ST]). 
Given a planar bipartite graph Q there are several Kasteleyn matrices K for Q but, as 
has been shown independently by David Wilson and Horst Sachs, the singular values of K 
or, equivalently, the eigenvalues of KK* , are independent of the choice of K. Following 
a question posed by James Propp ([P]), we search for a combinatorial interpretation for 
these numbers. 

In section 1 we introduce generalized Kasteleyn matrices for planar bipartite graphs 
and present a combinatorial interpretation for the determinant of such matrices A in terms 
of counting matchings. In section 2 we address the main issue of understanding what the 
coefficients of the characteristic polynomial of AA* represent, and then, in section 3, we 
consider the special case of Kasteleyn matrices. In section 4 we present the g-analogs of 
these ideas. In section 5 we take a look at rectangles in the plane. Finally, in section 6 we 
present a few other small examples. We find the language of homology theory helpful and 
use it throughout the paper. We thank James Propp, Richard Kenyon and Horst Sachs 
for helpful conversations and emails. 

1. Generalized Kasteleyn matrices and their determinants 

Let Q be a planar bipartite graph with n white vertices and n' black vertices. We 
number the white (black) vertices 1, 2, . . . , n (1, 2, . . . , n'). A generalized Kasteleyn matrix 
for Q is an n x n' complex matrix A such that 



I a i j I 



1, if the z-th white vertex and the j-th black vertex are adjacent, 
otherwise. 



Such matrices are conveniently represented by labeling the edges of Q with complex 
numbers of norm 1. 
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We may identify a generalized Kasteleyn matrix A with a cocomplex A e C 1 (^, S 1 ) 
by making the convention that indicates the value of A(ey ), being the oriented edge 
going from the j-th black to the i-th white vertex. A notational confusion must be avoided 
here: the complex numbers of norm 1 form a multiplicative group but the coefficients 
for homology or cohomology should be additive groups. Thus, from now on, the symbol 
S 1 shall denote the additive group R/Z and we denote the exponential x \— > exp(27rix) 
by rj : 8 1 — > C. In particular, we write = rj(A(eij)). Since C 2 (^,§ 1 ) = 0, any 
cocomplex A is automatically closed and a generalized Kasteleyn matrix A defines an 
element a E H 1 (^,§ 1 ). 

There is a natural inclusion Z/(2) C S ; this defines 77 : Z/(2) — > C with 77(771) = 
(-l) m . We also obtain induced inclusions C 1 (G,Z/(2)) C C 1 ^,^ 1 ) and H X {G, Z/(2)) C 
S 1 ). For a generalized Kasteleyn matrix A, A e C 1 (^,Z/(2)) if and only if A is a 
real matrix. 

Lemma 1.1: There is a unique element k e H X {Q, Z/ (2)) such that for any cycle C, 

k(C) = m + Z + l (mod 2) (1) 

where m is the number of vertices in the interior of C and 21 is the length of C. 

In the statement above, the word 'cycle' is used in the sense of graph theory: C is a 
simple closed curve in the plane composed of edges and vertices of Q and the interior of C 
is well defined by Jordan's theorem. Of course, graph theory cycles define homology cycles 
(i.e., closed elements of Ci(Q, Z)) but the converse is not always true; any homology cycle 
may nevertheless be written a linear combination of graph theory cycles. 

Proof: Uniqueness is obvious since the above equation gives the value of k computed 
against any cycle and thus, by linearity, against any element of H\(Q,Z) (recall that 
H \Q, Z/(2)) = Hom(ifi(0, Z), Z/(2))). 

In order to construct k, we first notice that H 1 ^, Z/(2)) = (Z/(2)) ?l , where h is 
the number of holes (bounded connected components of the complement) of Q: in this 
identification, the coordinates of a corresponding to a given hole is a(C), C being the 
outer boundary of the said hole. It is then clear that there exists a unique element of 
H 1 (Q 7 Z/(2)), which we call k, satisfying equation (1) for all such C. 

Figure 1.1 illustrates a minor complication which has to be kept in mind: the boundary 
of a hole is not always a cycle in the sense of graph theory. There should be no confusion, 
however: in (a) we have I = 2, m = 1 for the only hole and in (b) we have / = 2, m = 
and / = 2, m = 4 for the smaller and bigger hole, respectively. 




(a) (b) (c) 

Figure 1.1 
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Let C be an arbitrary cycle: we prove that equation (1) holds for C. The interior of 
C minus Q is a union of holes. If we discard the holes which are completely surrounded 
by other holes in C and consider the outer boundaries Ci, . . . , of the remaining ones, 
we have k(C) = ^^(Ci). Since equation (1) holds for each Ci we have 

k(C) = k + + Y mi (mod 2), 

i i 

where U and m>i correspond to Ci. Notice that L = I + Zj, L denoting the number of 
edges of some Ci on C or in the interior of C, and M = 21 + m — m^, M denoting the 
number of vertices of some Ci on C or in the interior of C. Finally, by Euler characteristic, 
k — L + M = 1 and we have equation (1) for C, proving our lemma. ■ 

We call k e H 1 ^, Z/ (2)) as defined in the previous lemma the Kasteleyn class; when 
the graph C? is not clear from the context, we write kg. The definition of k involves m 
and thus appears to depend on the way Q is drawn in the plane. Indeed, examples (a) and 
(c) in Figure 1.1 represent equivalent graphs, but the Kasteleyn classes are different; solid 
lines stand for a label 1 and dashed lines stand for a label —1. A Kasteleyn matrix is a 
generalized Kasteleyn matrix corresponding to the Kasteleyn class. 



We restrict ourselves for the rest of this section to the case n = n' in order to explore 
the relationship between matchings and the determinant of square generalized Kasteleyn 
matrices A. 

It is natural to interpret a matching of Q as the sum of its edges oriented from black 
to white and thus as an element of Ci(G, Z). The boundary of any matching always equals 
the sum of all white vertices minus the sum of all black vertices; thus, the difference of 
two matchings of Q is closed and may be identified with an element of Hi(Q,Z). Notice 
furthermore that the difference of two matchings may be written in a unique way as a sum 
of disjoint graph theory cycles. 

For a G H l (Q, S 1 ) and two matchings \i\ and [i2, v( a (P-2 ~ A*i)) is a complex number 
of absolute value 1. In particular, if a e H X (Q, Z/(2)) then r/(a(/i2 — A*i)) is 1 or —1. We 
then say that \i\ and \i2 have the same a-parity if r/(a(/X2 — A*i)) = 1; a-parity splits the 
set of matchings into two equivalence classes (occasionally one of these classes may turn 
out to be empty). 

Lemma 1.2: For any planar graph Q, for any a e H l (Q, S 1 ) and for any given matching 
Ho we have 



V{a(P2 - Hi)) = 



where Hi, A*2 and H range over all matchings. 
Proof: We may write the right hand side as 



5^ ?7(a(^ 2 -Ho))) [Yl V(^i ~ A*o)) J = [Yl ^( a (^2 - Ho))) \ Y ^( a (^o - Hi)) 



M2 



pi 



Ml 
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and distribute to get the left hand side. ■ 
Define 

s(&,g) = v(*(ij>i -^2)), 

where [i\ and ^2 range over all matchings; if Q admits no matchings we define S(a., Q) = 0. 
As an example, 5(0, Q) is the square of the number of matchings of Q. Also, for a e 
H X (Q, Z/(2)), 5(a, £/) is the square of the difference between the number of matchings in 
each a-parity equivalence class. 

A matching may also be though of as a bijection from the set of white vertices to the 
set of black vertices. Thus, if \i\ and \i 2 are matchings then \i^ x o /j 2 is a permutation 
of the set of white vertices: we say that these two matchings have the same permutation 
parity if and only if the this permutation is even. 

Lemma 1.3: Two matchings have the same permutation parity if and only if they have 
the same k-parifry, k being the Kasteleyn class. 

Proof: Let \i\ and \x 2 be two matchings and write (i\ — as a sum of disjoint cycles 
C\, . . . ,Cn of lengths 2Zi , . . . , 21n. The interior and exterior of any of these cycles is 
matchable, thus mi, . . . , tun as in Lemma 1.1 are all even. From equation (1), k(Cj) = h+l 
(mod 2) and thus k( / uo — £*i) = + 1) (mod 2). 

The permutation /j,q 1 o ^1 can be written as a product of N cycles (in the 
permutation sense) corresponding to Ci, . . . ,Cn with lengths Ii,...,In and the parity 
of the permutation ^ 1 o /j, 2 is thus J2(k + 1). This proves our claim. ■ 

Notice that permutation parity, unlike the Kasteleyn class, does not depend on how 
Q is drawn in the plane. A corollary of the previous lemma is thus that if differences 
of matchings generate Hi(Q,Z) then the Kasteleyn class of Q is the same for all planar 
embeddings. 

Lemma 1.4: For any generalized Kasteleyn matrix A we have \ det(A)| 2 = <5(a + k, Q). 

Proof: Each non-zero monomial in the expansion of det(A) corresponds to a matching. 
Thus, each matching \i contributes with a complex number of absolute value 1 to det(A). 
The expression n(a(/i — /i )) obtains, up to a fixed multiplicative constant of absolute 
value 1, the product of the corresponding elements of A. From Lemma 1.2, k-parity is 
permutation parity, i.e., gives the sign of the monomial in the definition of the determinant. 
Thus, the contribution of (i to det(A) is, again up to a fixed multiplicative constant of 
absolute value 1, rj((a + k)(p — no)), proving our lemma. ■ 

As a special case, if K is a Kasteleyn matrix, | det(_ftT)| is the number of matchings of 
Q: this is Kasteleyn's original motivation. 

2. Singular polynomials of generalized Kasteleyn matrices 

Having provided an interpretation for | det(A)| when A is square, it is natural to ask 
about other functions of A, specially if A is not square. We should not expect natural 
interpretations for the argument of det(A) since it depends on the way we assign labels 
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to vertices. Also, a few simple experiments will show that the spectrum of A (even if A 
is square) is not a function of a. The following lemma tells us what functions of A are 
determined by a. 

Lemma 2.1: Let A be a generalized Kasteleyn matrix for Q and let a be the corresponding 
element of i? 1 (^,S 1 ). Then the generalized Kasteleyn matrices for Q also corresponding 
to a are precisely the matrices of the form DiAD 2 where Di and D 2 are unitary diagonal 
matrices. Furthermore, if Q is connected, D\AD 2 = D' 1 AD' 2 if and only if there exists a 
complex number z of absolute value 1 with D\ = zD[, D 2 = z~ x D' 2 . 

It is possible to give a more elementary proof, but following the spirit of the rest of 
this paper we phrase the proof in homological language. 

Proof: As we saw in Section 1, generalized Kasteleyn matrices correspond to 1- 
cocomplexes in C 1 (^,S 1 ); two such 1-cocomplexes A and A' induce the same element of 
H X {Q, S 1 ) if and only if their difference is the coboundary of a 0-cocomplex. A 0-cocomplex 
D is a function assigning an element of S 1 to each vertex; the ry's of these elements may 
conveniently be arranged in a pair of unitary diagonal matrices, D w for the white and 
Dt, for the black vertices. It is a simple translating process to verify that the cocomplex 
A + <i(D) corresponds to the generalized Kasteleyn matrix D W AD^ X , thus proving our first 
claim. The uniqueness of Di and D 2 up to a constant multiplicative factor corresponds to 
the fact that the only closed O-cocomplexes are the constants, i.e., that rj"°(^,S 1 ) = S 1 if 
Q is connected. ■ 

We recall that for any complex nx n' matrix B, there are unitary matrices Ui and U 2 
such that S = U\BU 2 is a real diagonal matrix with non-increasing non-negative diagonal 
entries; S is well-defined given B and its diagonal entries (i.e., the su entries of S, even if 
S is not square) are called the singular values of B. The rows of U\ (resp., columns of U 2 ) 
are called the left (resp., right) singular vectors of B. It is easy to see that the singular 
values and left (resp. right) singular vectors of B are the non-negative square roots of the 
eigenvalues and the eigenvectors of BB* (resp., B*B). Inspired in these classical notions, 
we call the characteristic polynomial of BB* the singular polynomial of B: its roots are 
the squares of the singular values of B. Also, the singular polynomials of B and U\BU 2 
are equal and the singular polynomials of B and B* differ by a factor of t n ~ n . 

It follows from Lemma 2.1 and the remarks in the previous paragraph that the singular 
polynomial of A is determined by a: we call it P a . The singular values of A and, if the 
singular values are simple, the absolute values of the coordinates of the singular vectors 
(up to a constant factor) are also determined by a. We shall now present what we find 
to be a reasonably natural interpretation for the coefficients of P a . While these numbers 
determine the singular values the question remains whether a nice interpretation exists for 
the actual singular values and vectors. 

Let 7i C Q be a balanced subgraph of Q: the inclusion induces a map ng^n : 
rJ 1 (^,S 1 ) — > H 1 (7i,S 1 ). More concretely, if a corresponds to a generalized Kasteleyn 
matrix A then Tig ^(a) corresponds to the submatrix of A obtained by picking only the 
elements for which both row and column correspond to elements of H. The simplest 
interpretation is probably in terms of labels for edges: just keep the old labels. When this 



5 



causes no confusion, we write a instead of TTg^a): for instance, we write S(a., Ti) instead 
of the more correct but cumbersome S(irg ^(a)? H). 

Theorem 2.2: Let A be a generalized Kasteleyn matrix and let P a (t) = t n + ait n_1 + 
■ • ■ + a n -it + a n be the singular polynomial of A. Then 

o m = (-l) m *(a + k w ,W) (2) 

\H\=2m 

where 7i ranges over all balanced subgraphs with 2m vertices. 

Notice that for m = n equation (2) is equivalent to Lemma 1.4. For m = 1, we get the 
simple remark that |ai| is the number of edges of Q, regardless of a. An interpretation for 
ci2 is already subtler: each subgraph with two white and two black vertices contributes with 
a real number between and 4. Subgraphs which are not matchable of course contribute 
with and two disjoint edges as well as four points on a line contribute with 1. The 
interesting part are the squares, which admit two matchings, say Hi and \i2- then 5(a + 
k-H, H) = |1 — t](a(ni — /x 2 ) ) | 2 ; in general, this may be any number between and 2 but if 
a G if 1 (^,Z/(2)) then this is or 4. 

In order to prove this Theorem, we need an auxiliary result in linear algebra. The 
proof of Lemma 2.3 (actually a rather straightforward computation) may be found in [??]. 

Lemma 2.3: Let P(t) = t n + ait n_1 H h a n -it + a n be the singular polynomial of A 

(where A is an arbitrary n x n' complex matrix). Then 

a m = (-irj2\ detB \ 2 

B 

where B ranges over all m x m submatrices of A. 

Proof of Theorem 2.2: Since balanced subgraphs of Q with 2m elements correspond 
to m x m submatrices of A, this is a consequence of Lemma 1.4 and Lemma 2.3. ■ 

We now describe another, more graphical, interpretation for Theorem 2.2. We define 
a pipe system of Q as an oriented pair v = (^i, ^2) of matchings of a subgraph H of Q; we 
call 7i the support of the pipe system. Figure 2.1 shows an example of a pipe system: we 
draw the edges of \ii oriented from black to white and the edges of [i\ from white to black 
(unused edges are represented by dotted lines) . A pipe system is thus a collection of pipes 
(i.e., oriented edges of Q) such that, at each vertex, there is either one pipe coming in and 
one pipe going out or no pipe coming in and no pipe going out (the water that comes in 
must go out and you can not pipe too much water through a vertex). We define the size 
\u\ of the pipe system as half the number of vertices in TL (in Figure 3.3, \v\ =5). 

The Kasteleyn class shall be called k^. A pipe system obtains an element of 
Ci(7i,Z) (and thus of Ci(C/,Z)) but must not be confused with it: if two pipes cancel 



6 



Figure 2.2 



each other home-logically, they still have to be taken into account for the pipe system. If 
a G H 1 (H, S 1 ), 7](&(is)) and //((a+ k„)(i/)) are well defined complex numbers. 

Corollary 2.4: Let A be a generalized Kasteleyn matrix and let P a (t) = t n + ait 71-1 + 
■ ■ ■ + a n -it + a n be the singular polynomial of A. Then 

^a(t) = ^(-l) IH t n - |H r 7 ((a + k ! ,)(z,)), 
where v ranges over all pipe systems of Q . 

Proof: This follows directly from Theorem 2.2 and the definitions. ■ 

3. Singular polynomials of planar graphs 

Theorem 2.2 provides an interpretation for the coefficients of singular polynomials of 
arbitrary generalized Kasteleyn matrices. In this section we take a closer look at the right 
hand side of equation (2) when A is a Kasteleyn matrix. 

For planar balanced bipartite graphs 7i C define pg,n € H 1 (H : Z/(2)) by Pg,n = 
7Tg ; -ft(kg) — kft. In Figure 3.1 we illustrate the several objects involved in this definition: 
(a), (b), (c) and (d) represent kg, Tvg^(kg), k^ and pg,n, respectively, where again solid 
lines stand for a label 1 and dashed lines stand for a label —1. The following lemma 
provides an alternate definition for this class. 

Lemma 3.1: Let H C Q be balanced planar graphs and let C be a cycle in H. Let q be 
the number of vertices of Q not belonging to H which are inside C. Then 

PG,n( C ) = Q ( mod 2 )- 



Proof: This follows directly from equation (1) in Lemma 1.1. ■ 




(a) (b) (c) (d) 

Figure 3.1 
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In the hope of making the intuitive meaning of this definition clearer, especially for 
adjacency graphs of quadriculated or triangulated disks, we introduce some extra structure. 

Let Q be the CW-complex obtained from Q by closing each hole with a 2-cell; Q is thus 
always homeomorphic to a disk. For Ti C let Ti C Q be the obtained from Ti by adding 
the 2-cells of Q whose boundaries are contained in TC; in other words, we close the holes 
of 7i which contain no points of Q. The inclusion 7i C H induces an injective map from 
H 1 (H,Z/(2)) to H 1 {H,1/{2)) which allows for a natural identification of H 1 (H,Z/(2)) 
with a subset of H l (H, Z/(2)). 

Lemma 3.2: p g>w belongs to H X {H, Z/(2)). 

Proof: This follows easily from Lemma 3.1. ■ 

Recall that two tilings by dominoes of a quadriculated region are said to differ by a Hip 
if they coincide except for two dominoes; in other words, their difference (in the homological 
sense) is a square. If Q is the graph of a quadriculated planar region, the difference between 
two tilings of Ti differing by a flip is in Hi(H,Z); thus, tilings mutually accessible by 
flips always have the same pg^-parity. 

We may now state the promised interpretation for the coefficients of singular 
polynomial P^ of K. Since P^ is well defined from Q, we may adopt a lighter notation and 
call it Pg, the singular polynomial of Q. 

Theorem 3.3: Let Q be a planar bipartite graph and let Pg = t n +kit n ~ 1 -\ \-k n -it+k n 

be the singular polynomial of Q. Then 

h m = (-l) m 6 (P^n,n) (3) 

\H\=2m 

where Ti ranges over all balanced subgraphs with 2m vertices. 

Proof: This is a corollary of Theorem 2.2 and the definition of pg t n- • 

Recall that if Pg,n = then d(pg t -n, Ti) is just the square of the number of matchings 
of TL. This always happens if TC is simply connected. As a corollary, if Q is the graph of a 
quadriculated planar region and m < 3, or if Q is the graph of a triangulated planar region 
and m < 5, then 

k m = (-ir Yl 5 ^ H ) 

\H\=2m 

where H ranges over all balanced subgraphs (subregions) with 2m vertices (squares, 
triangles). 

Notice that pg^ an d thus the right hand side of equation (3), depends on the 
way Q is drawn in the plane. Examples (a) and (c) in Figure 1.1 show that Pg indeed 
depends on the way Q is drawn: for (a) we have hi = 9 but for (b) we have hi = 5. 
This causes the singular values to change in a complicated way: for (a) the singular 
values are approximately 0.5549581321,0.8019377358,2.246979604 while for (c) they are 
0.3472963553, 1.532088886, 1.879385242. 
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It is natural to conjecture that the number of non-zero singular values coincides with 
the size of a maximal partial matching of Q. In Figure 3.2(a) we present an example 
to show that this is not always true: there are partial matchings of size 3 but since the 
singular polynomial of Q is t 3 — 7t 2 + lOt there are only two non-zero singular values: a/2 
and a/5. In Figure 3.2(b) we draw the same graph in a different way and we now have 
three non-zero singular values: 1, \/2 and 2. 




(a) (b) 
Figure 3.2 



We state Theorem 3.3 in the language of pipe systems. Denote pg^i^) (where 7i is 
the support of v) by p(^) e 2/(2). We describe an elementary definition of p{y). J°i n 
pairs of vertices not in Ti, matching black vertices with white vertices in an arbitrary way; 
if all intersections are transversal, p(f) is the parity of the number of intersections of such 
new lines with the pipes. 

Corollary 3.4: Let Q be a planar graph and Pg(t) its singular polynomial. Then 

Pg(t)=J2(- 1 ) Htn ~ H v(p(")), 
where v ranges over all pipe systems of Q . 

Proof: This follows directly from Theorem 3.3 and Corollary 2.4. ■ 
4. q 

In several branches of combinatorics, (/-analogues or quantizations of classical problems 
have been seen to be interesting and useful. There are often several interpretations for the 
(/-analogue of a given concept, some sophisticated (involving quantum groups and the 
like) and some elementary. In this section we briefly consider a (/-analogue of Kasteleyn 
matrices in a very naive way and extend the results of the previous sections to this setting; 
our interest in doing so is that the methods of the previous sections extend very easily to 
this more general context and the coefficients will actually have a natural interpretation. 

Let Cq = C[(/,(/ -1 ]. We extend the usual complex conjugation to C q by postulating 
q = (/ _1 ; q may be thought of as an unknown complex number of absolute value 1. Let 
Q be a planar bipartite graph with n white vertices and n' black vertices. A generalized 
Kasteleyn q-matrix for Q is an n x n! matrix A with coefficients in C q such that 

J 1, if the i-th white vertex and the j-th black vertex are adjacent, 
a ij a tj - | q otherwise; 
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thus, the entries are always monomials and substituting q by a complex number of absolute 
value 1 changes a generalized Kasteleyn (/-matrix into an ordinary generalized Kasteleyn 
matrix. 

Consider the additive group § g = S 1 © Z and let q be the canonical generator of the 
Z component. If we extend the classical rj to rj : S 1 © Z — > C q by postulating 77 (q) = g we 
may identify a generalized Kasteleyn (/-matrix A with a cocomplex A e C X {Q, S 9 ). Again, 
C (£, §,) = and A defines an element a e fl" 1 ^, §,). 

Let Q be a planar graph. Bounded connected components of the complement of Q have 
well defined positively oriented boundaries (3 in H\(Q, Z). We define a Kasteleyn (/-matrix 
to be a generalized Kasteleyn (/-matrix A such that a(/3) = q for all such boundaries /3. We 
define the singular q-polynomial of £ to be the singular polynomial of a Kasteleyn (/-matrix 
of Q. As before, singular (/-polynomials are easuly seen to be well defined but now they 
are of course polynomials in C q [t], or, rather equivalently, polynomials in two variables 
q and t. Finally, define the area of a pipe system u, A{y) to be the number of bounded 
connected components of the complement of Q positively surrounded by is, counted with 
sign and multiplicity. 

With these definitions we have the following theorem. 

Theorem 4.1: Let Q be a planar graph and Pg(q, t) its singular q-polynomial. Then 

P 6 (q,t) = J2(- 1 ) H <l A{u) t n ~ H v(p^)), 

where v ranges over all pipe systems of Q . 

Since the proof is entirely analogous to that of Corollary 3.4, we leave the details to 
the reader. 



5. Rectangles 

Kasteleyn ([K]) computes the determinant of K for rectangles essentially by computing 
its singular values. For the reader's convenience, we repeat that part of his work in our 
language. In order to simplify notation in the statement and proof, let 

X+^ N = e I 1 < k < an d if k = then 1 < i < 

X~ N = \ (k, i) E Z 2 I 1 < k < M ~ H 1 and if k = M ^ 1 then 1 < £ < 



Theorem 5.1: Let Q be a M x N rectangular grid and let K be its Kasteleyn matrix. 
Then the non-zero singular values of K are (7k,t, (k, £) € X^ N , where 

al t = {a k + a- k Y + ^ + r*)\ a = exp , (1 = exp . 
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The complicated description of the allowed values of the indices k and £ is necessary 
in order to avoid zeroes and duplications in a way which is correct for all possible parities 
of M and N (Kasteleyn has a simpler formula since he assumes N to be even). Notice 
that 



kn 9 £n 



Oh o = 2 \ cos 2 V cos 2 



M + l N+l 



1/2 



(an expression closer to Kasteleyn's), a kji = a M +i-k,t = cr k ,N+i-i = cr M +i-k,N+i-i and 
that a k ,e = if and only if M and N are both odd, k = (M + l)/2 and £ = (N + l)/2. 
Thus, if we just demand 1 < k < M and 1 < £ < N then all non-zero singular values are 
counted twice and we occasionally introduce a 0. 

Proof: We index vertices by pairs (k',£'), 1 < k' < m, 1 < £' < n. The vertex (k',£') 
is called white when k' + £' is even. Define K as the Kasteleyn matrix with entries 1 for 
horizontal edges and i for vertical edges: K defines a linear transformation from the "black 
space" to the "white space" . Consider the white vectors 

Wk,i — (« - a ){P - P ), 

(k 7 £) E X+y. they clearly form an orthogonal basis for the white space (this is where a 
careful choice of X+^ n becomes necessary). Similarly, the black vectors bk,i defined by the 
same formula with (k,£) G n form an orthogonal basis for the black space. A simple 
computation yields \wk,i\ = \bk,l\ f° r G X^^ n and 

Kb k , t = {{a k + a~ k ) + i((3 e + p~ £ )) w k , e , 
K*w k/ = ((a k + a~ k ) - i(p £ + (3~ e )) h, e . 

Thus, w kj £ and b k ,e are singular vectors and a k ^ are singular values. ■ 
Corollary 5.2: Let Q be a m x n rectangular grid. Let 



a = ex P ( — ~r~r ) > /3 = exp ( — ^— ) , N 



mn 



Then 



J! (t - (a k + a~ k f - ((3* + (3~ e ) 2 ) = E ^(-1)* ^ *(P*,w,W), 

(k,t)€X~ in 3=0-N \H\=2j 

where 7i ranges over all balanced subgraphs with 2j vertices. 

Proof: This follows directly from Theorem 3.3 and Theorem 4.1. ■ 

These results show that the characteristic polynomial of KK* usually factors a lot if 
Q is a rectangle. If £ is a root of unity whose order M is the least common multiple of 
2(m + 1) and 2(n + 1) then all the roots cr 2 e of this polynomial are in M D Z[£], a ring of 
degree <f>(M)/2 over Z. In particular, for square grids of order n, irreducible factors of the 
characteristic polynomial of KK* have degree at most n. Here are a few sample examples; 
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we give the polynomial det(tl — KK*) = t n + kit n 1 + h /c n _it + k n (whose roots are 

the squares of singular values) factored in Z. 



[1,1] 1 

[2, 1] t-1 

[2,2] (t-2) 2 

[3, 1] t-2 

[3,2] (t-1) (t-3) 2 

[3,3] (t-2) 2 (t-4) 2 

[4,1] t 2 — 3t + 1 

[4,2] {t 2 -ht + h) 2 

[4,3] (t 2 -3t + l) (t 2 -7t + ll) 2 

[4,4] ( t 2_ 6t + 4 ) 2 (t _ 3) 4 
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{t-l){t-2f{t-Af 

(t -l)(t- 2) (t - 5) 2 (* - 3) 3 

(t 2 -3t + l) (t 2 - 5t + 5) 2 (t 2 - 9t + 19) 2 

(t - l) 2 (t - 2) 2 (t - 3) 2 (* - 6) 2 (t - 4) 4 

t 3 -5t 2 + 6t-l 

(t 3 -8t 2 + 19t-13) 2 

(t 3 -5t 2 + 6t-l) (t 3 - lit 2 + 38t - 4l) 2 

(t 6 - 19t 5 + 142t 4 - 529t 3 + 1017t 2 - 922t + 28l) 2 

5t 2 + 6t - l) (t 3 - 8t 2 + 19* - 13) 2 (t 3 - Ut 2 + 63t 

(t 3 - 10t 2 + 24 1 - 8) 2 (t 3 - 10t 2 + 31 1 - 29) 4 

13 



[7,7] 



(t - 2) 2 (t 2 - At + 2) 2 (t 2 - 8t + 8) 2 (t 2 - 8t + 14) 4 (t - 4) 



[8, 8] (t - 2) 2 (t 3 - 12t 2 + 36t - 8) 2 (t 3 - 9t 2 + 2At - 17) 4 (t 3 - I2t 2 + 45t - 53) 4 
6. Examples 

Although Aztec diamonds have so many interesting properties (see [EKLP] and [P]), 
the characteristic polynomial of KK* does not factor very much: 



2- Aztec diamond (t s - 8t 2 + 17t - 8) 



3-Aztec diamond (t 4 - 10t 3 + 28t 2 - 2At + A) 2 (t - 4) 4 



4- Aztec diamond 

(t 10 -32t 9 +441t 8 -3424t 7 +16432t 6 -50240t 5 +97041t 4 -112896t 3 +70921t 2 -18784t+1024) 2 



(t 11 - 34t 10 + 496t 9 - 4064t 8 + 20562t 7 - 66524t 6 + 137728t 5 
5- Aztec diamond , „ n 

-177120t 4 + 131825t 3 - 49066t 2 + 6576t - 128) 2 (t - 4) 8 

The fact that these polynomials are always squares follows from symmetry. The 
factor (t — 4) 4k seems to appear in the 2k + 1- Aztec diamond, a fact for which we have no 
explanation. 

Finally, here is a small list of "irregular" examples. A possible real Kasteleyn matrix 
is indicated by the dashed lines (the — l's). 



t 3 - 7t 2 + 15t-9 = (t- l)(t - 3) 2 
1.732051,1.732051,1 



t 3 -m 2 + m-4 = (t-4)(t-i) 2 

2,1,1 
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t 4 - 10t 3 + 35t 2 - 50t + 25 = [f -5t + 5) 2 
1.902113, 1.902113, 1.175571, 1.175571 



t 4 - 8t 3 + 20t 2 - 16t + 4 = (t 2 -4t + 2) 2 
1.847759, 1.847759, 0.765367, 0.765367 




t 5 - 13t 4 + 63t 3 - 140t 2 + 140* - 49 = 
(t 2 - 6t + 7)(t 3 - 7t 2 + Ut - 7) 
2.101003, 1.949856, 1.563663, 1.259280, 0.867768 




t 5 - 13t 4 + 62t 3 - 132t 2 + 121* - 36 = 

(t - 4) (t 4 - 9t 3 + 26t 2 - 28t + 9) 
2.126757, 2, 1.576415, 1.197126, 0.747468 
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